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We discuss new approaches to compute numerically the Casimir interaction 
energy for waveguides of arbitrary section, based on the boundary methods 
traditionally used to compute eigenvalues of the 2D Helmholtz equation. These 
methods are combined with the Cauchy's theorem in order to perform the 
sum over modes. As an illustration, we describe a point-matching technique 
to compute the vacuum energy for waveguides containing media with different 
permittivities. We present explicit numerical evaluations for perfect conducting 
surfaces in the case of concentric corrugated cylinders and a circular cylinder 
inside an elliptic one. 
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1. Introduction 

In this paper we will be concerned with the numerical calculation of the 
Casimir interaction energy in geometries with translational invariance along 
one direction, i.e. very long cylinders of arbitrary section. For the sake of 
simplicity, we will first discuss the case of a massless quantum scalar field 
that satisfies Dirichlet or Neumann boundary conditions on the surfaces of 
the cylinders. As we will see, in some particular situations the generalization 
to the electromagnetic field and/or more general boundary conditions will 
be straightforward. 

The Casimir energy is formally given by 

EM = lim \ Y.ie-'-'-^Wp - e-'^^^w,) , (1) 

P 

where Wp are the eigenfrequencies of the scalar field satisfying the appro- 
priate boundary conditions on the surfaces of the shells, and Wp are those 
corresponding to a situation in which the distances between the shells is 
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very large. The subindex p denotes the set of quantum numbers associated 
to each eigenfrequency. We have introduced an exponential cutoff for high 
frequency modes. 

For the particular geometry considered here, p = [n, k^) and the eigen- 
frequencies are of the form uin,k^ = -^/fc^ + A^, where kz is a continuous 
variable associated to the translational invariance along the z-direction and 

are the eigenvalues of the Laplacian on the two-dimensional transversal 
section E contained in the plane (x,y): 

A^u = -\^u . (2) 

The eigenfunctions u(x) satisfy Dirichlet or Neumann boundary conditions 
on r, the boundary of S. 

The Helmholtz equation (2) arises in many branches of physics, from 

the vibration of membranes to quantum billiards, and there arc a plethora 
of methods to compute numerically its eigenfunctions and eigenvalues.^ 
Among them, the "boundary methods" are based on the following strategy: 
the solution u is written as a (finite) linear combination of basis functions 
that satisfy Helmholtz equation inside S. The coefficients of the linear com- 
bination are chosen in such a way that the boundary conditions are satisfied 
at a finite number of points on F. The linear system of equations that de- 
termine the coefficients has a non trivial solution only for some particular 
values of A, the eigenvalues of the system. 

For example, in the Point Matching Method (PMM),^ one expands the 
eigenfunction u in terms of a basis of solutions of the Helmholtz equation 
in free space (p^^^ (x) 

«(x)=f;a,^f)(x). (3) 
j=i 

In the numerical calculation this expansion is truncated at given j = N, and 
the boundary conditions are imposed on A'^ points on F. These boundary 
conditions become a set of homogeneous, linear equations for the unknown 

coefficients aj (Ma = 0, with M a A-dependent N x N matrix) which has 
nontrivial solutions only when detM = 0. The last equation can be used to 
determine numerically the eigenvalues A„. 

In a similar approach, known as the Method of Fundamental Solu- 
tions (MFS),'^ the eigenfunction u is expanded in terms of solutions of the 
Helmholtz equation with a point source at an arbitrary location Sj , that we 
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denote by ux{x, sj) 

CO 

If the sources are located outside S, this is a solution of the homogeneous 
Helmholtz equation inside S. Once again, the sum is truncated at j = N, 
and the coefScients bj are determined by solving the linear system that re- 
sults when imposing the boundary conditions on a finite number of points 
on r. The roots of the determinant of the associated matrix are the eigenfre- 
quencies of the problem. This is the simplest version of the MFS, in which 
the locations of the sources are fixed. 

One can find in the literature discussions about spurious solutions, im- 
provements and alternative methods to find the eigenvalues. We refer the 
reader to Refs.^'^ for more details. 

A crucial point is that, at a practical level, the knowledge of the spec- 
trum of the Helmholtz equation is not enough to compute the Casimir 
energy. The reason is that the numerical evaluation of the sum over modes 
in Eq.(l) is extremely unstable,^ and one has to subtract very large num- 
bers to compute the finite interaction energy. The calculation is complicated 
even for the simplest case of Casimir effect in 1 + 1 dimensions. 

Instead of performing explicitly the summation, it is far more efficient 
to combine the methods mentioned previously with the Cauchy's theorem 

where f{z) is an analytic function in the complex z plane within the closed 
contour C, with simple zeros at zi,Z2,--- within C. We use this result 
to replace the sum over the eigenvalues of the Helmholtz equation in the 
Casimir energy Eq.(l) by a contour integral. In this way, it is not necessary 
to solve numerically the equation detM = for the eigenvalues, but to take 
/ = detM in the Cauchy's theorem. In other words, if in the numerical 
method to solve the Helmholtz equation the eigenvalues are the roots of 
a given function, one can integrate this function in the complex plane in 
order to get the Casimir energy. The combination of the use of numerical 
methods to compute the eigenvalues with the Cauchy's theorem is the main 
idea we want to put forward in this paper. 

In the next Section we will describe the simplest version of the PMM to 
a situation in which the surfaces separate regions of different permittivities, 
generalizing our previous results^ for perfect conductors. In Sections 3 and 4 
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we will review some numerical evaluations of the Casimir interaction energy 
for perfect conductors. In Section 5 we include our final remarks. 

2. Point-Matching Numerical Approach 

A media-separated waveguide presents an interesting setup for the applica- 
tion of the PMM. This technique has been widely used to solve eigenvalue 
problems in many areas of engineering science.^ The boundary conditions 
are imposed at a finite number of points around the periphery of both 
media. 




Fig. 1. A two-separated media waveguide in which one conductor encloses two diflferent 
dielectric media. Each has arbitrary cross section. 

For the sake of concreteness, we will bear in mind the situation in which 
one perfect conductor encloses two dielectric media, as shown in Fig.l, 
although the method could be applied to more general cases. 

The general solution of the Helmholtz equation in region I (inside the 
inner cylinder) is 

w = ^^„J„(A(^V)e*'»^ (6) 

m 

while in region II (annular region) 

u = Y^[BmJm{X^"^r) + CmHj^\X^''\)]e'"'', (7) 

where (r, 0) are polar coordinates, and Jm and Hm^ are the m-th order 
Bessel functions. The constants Am,Bm and Cm are determined by the 
boundary conditions. In both equations we have defined A*^''-' = ta'^'^ — k^. 
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We assume the outer surface to be a "perfect conductor", and impose 
Dirichlet boundary conditions on a finite number of points (r,, 9q) of C2: 
s 

0= [B„J„(A(^^)r,)+C„ifW(A(^^)r,)]e^'"^', (8) 

m=-S 

(alternatively, for the TE modes of the electromagnetic field, one should 
impose Neumann boundary conditions). The surface Ci as a dielectric in- 
terphase separating media ei and €2, and therefore we impose continuity of 
the field and its derivative: 
S S 
^ A„J„(A(^)rp)e'™^- = 5^ [S„J„(A(^^)rp)+C„if«(A(")rp)]e'™^- (9) 

m=—S m——S 
S UII) s 

J2 AMx(%,)e^^^'' = ^ J2 [S™J;(A(^'V,) + ^7„^f;(^)(A(")r,)]e*™^^ 

m——S m——S 

where {rp,6p) are points on the curve Ci. 

The boundary conditions can be written, in matrix form, as 

= NiB + N2C, 
RiA = MiB + M2C, 

R2A = M[B + M'^C. (10) 

Eliminating the coefficients we end with 

NiB + N2C = Q, PiB + P2C = 0, (11) 

where Pi and P2 can be written as 

Pi=Mi-RiR^^M[, P2 = M2-RiR2^M^. (12) 

It is worthy to note that as RiRl^^ is proportional to A2/A1, then 
i?ii?2^^ — > when ei — > 00. Thus, the matrices Pi — > Mi and P2 — >■ M2, re- 
obtaining in this way, the usual perfect conductor wave-guide case studied 
in.6 

For the system of Eq.(ll) to have non trivial solutions, the determinant 
must be zero, i.e. 



det 



Ni N2 
Pi P2 



= detP2 ■ detNi ■ det{l - iV2P2"^-Pi^r^) = . (13) 



This equation determines the eigenfrequencies associated to the geometry. 
However, as already mentioned, in order to compute the Casimir energy it 
is not necessary to find each eigenvalue but to integrate the determinant 
Q = det{l — N2P2^PiNi^) in the complex plane. 
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We have developed a numerical Fortran routine in order to evaluate the 
Casimir interaction energy in the case in which the field satisfies Dirichlet 
or Neumann boundary conditions on both curves Ci and C2. In this case 
one should consider the fields only in region II with £2 = 1- After some 
straightforward steps one can re-write the Casimir energy as a single in- 
tegral in the imaginary axis iy = A^^^^ . For Dirichlet boundary conditions 
the result is 

L f°° 

^12 = ^ dyylnQiiy), (14) 

while for Neumann boundary conditions one can derive a similar expression 
with a different function Q. It is worth to stress that these Casimir energies 
correspond to those of TM and TE modes of the electromagnetic field in 
the presence of perfect conductors. 



3. Cylindrical rack and pinion 

When two concentric cylinders have corrugations, the vacuum energy pro- 
duces a torque that could, in principle, make one cylinder rotate with re- 
spect to the other. This "cylindrical rack and pinion" has been proposed in 
Ref.,'' where the torque has been computed using the proximity force ap- 
proximation. It was further analyzed in,^ where the authors obtained per- 
turbative results for Dirichlet boundary conditions in the limit of small am- 
plitude corrugations. In this Section, we numerically evaluate the Casimir 
interaction energy for two concentric corrugated, perfect conductor cylin- 
ders. The cylinders have radii a and b, and we will denote by r_ = b — a 
the mean distance between them and by r+ = a + b the sum of the radii. 
We will use the notation a = b/a. The points in the mesh, that give us the 
corrugated cylinder boundaries, are described by the following functions: 

ha{0) = h sm{iy0) ; hbi9) = h sm{i^9 + 0o), (15) 

where h is the corrugation amplitude and u is the frequency associated 
with these corrugations. The Casimir torque can be calculated by taking 
the derivative of the interaction energy with respect to the shifted angle 

In Fig. 2 we show the numerical evaluation of the TM Casimir interaction 
energy for this geometry. The plot shows the results obtained using the 
PMM with a = 2 and corrugation frequency i' = 3, for different values 
of the amplitude of the corrugation h. As expected the amplitude of the 
oscillations grows with h. For each value of h we have performed a numerical 
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fit of the data in order to compare with the analytical prediction. With 
dotted lines we have plotted the fit y{x) = A * cos{x) for each curve in 
Fig. 2. The agreement between dots and dotted lines is extremely good. 
Similar results can be obtained for the Neumann (TE) modes (see® for 
details). 
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Fig. 2. Casimir interaction energy (TM modes) as a function of <pQ for a = 2 and 
different values of the perturbation h. The different shaped dots are the numerical data 
obtained with our program while the lines represent the numerical fit of each curve. 
Energies are measured in units of L/a^, and distances in units of a. 



It is worth to remark that, when the amplitude of the corrugation is 
not very small, the exact results cannot be reproduced with a simple fit of 
the form y{x) = A* cos(a;). This is illustrated in Fig. 3, where we see that, 
for the biggest corrugated amplitude h = h/a = 0.3, the exact result differs 
from the cosine function.® 



4. Outer conductors with focal lines: cylinder inside an 
ellipse 

Some time ago, there was a conjeture^ based on a geometric optics approx- 
imation, about the possibility of focusing vacuum fiuctuations in parabolic 
mirrors. It was argued that a parabolic mirror is capable of focusing the 
vacuum modes of the quantized electromagnetic field, therefore creating 
large physical effects near the mirror's focus. With this motivation, in this 
Section we shall evaluate the Casimir interaction energy for configurations 
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Fig. 3. Casimir interaction energy (TE and TM modes) as a function of (^o for a = 2, 
1/ = 3 and h = 0.3. The different shaped dots are the numerical data obtained by our 
program while the line represents the numerical fit of each curve. In this case, the plot 
shows that the exact result cannot be fitted by a function y{x) = A* cos(a;). Energies 
are measured in units of L/a?. 



in which the outer conducting shell has a cross section that contains focal 
points. 

We will consider one small inner cylinder and an outer ellipse. We will 
denote by a the radius of the inner cylinder, by bi and 62 the minor and 
major semiaxes of the ellipse, respectively, and by / the distance between 
the foci and the center of the ellipse. The coordinates of the center of the 
cylinder with respect to the center of the ellipse will be {ex, Cy)- We will use 
an additional tilde to denote adimensional quantities, i.e distances in units 
of a: h = bi/a,f^f/a, etc. 

For this configuration, we use a mesh where with an inner cylinder, and 
an outer ellipse with semiaxes hi — 4 and 62 = 4.33. The ellipse has two 
focal points at / = 1.66. We present the results for the Casimir energy in 
Fig.4. 

From Fig.4 it is possible to see that there is an unstable equilibrium 
position at the origin under displacements of the inner cylinder along the 
(vertical) Ey direction. As expected, it is also possible to check that the 
energy grows as well as the cylinder gets closer to the surface of the outer 
ellipse. Fig.4 also shows a monotonic behaviour of the energy as a function 
of the position, even when passing through the focus. So we do not see a 
focusing of vacuum fluctuations near the focus of the ellipse. However, in 
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Fig. 4. Numerical evaluation of the Casimir interaction energy for an inner cylinder an 
eccentric outer ellipse, as a function of the position of the cylinder along the vertical 
axis. Energies are measured in units of L/o? . 



order to confirm this result one sliould consider much smaller inner cylin- 
ders, in order to explore shorter wavelengths. This will require much more 
computational effort. 

Finally, we have also checked^ that there is an unstable equilibrium 
position at the origin when moving the inner cylinder in the (horizontal) 
err direction. 



5. Final remarks 

We have presented new numerical methods to compute the vacuum en- 
ergy for arbitrary geometries with translational invariance. The approach 
is based on the use of traditional boundary methods to compute eigenval- 
ues of the two dimensional Helmholtz equation, combined with Cauchy's 
theorem. 

As a particular example, we have described a straightforward version 
of the point-matching method to compute the Casimir interaction energy 
for a waveguide with different permittivities, and reviewed some numerical 
calculations for perfect conductors. In all examples, for the numerical cal- 
culations we have chosen pair of points with the same angular coordinate 
with respect to the inner cylinder. For less symmetric configurations, and 
when the surfaces of both conductors are closer to each other, it will be 
necessary to consider grids with a larger number of points, and to optimize 
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their positions. As in the appUcations to acoustic or classical electromag- 
netism, special care must be taken for surfaces with pronounced edges, 
clefts or "handles", where the point-matching technique may not be accu- 
rate to determine the cigcnfrequencics. In these cases, more sophisticated 
approaches'' could be necessary to optimize the numerical evaluation and 
to avoid spurious solutions. 

We would like to thank Kim Milton for the organization and his kind 
hospitality during QFEXT09. This work has been supported by CONICET, 
UBA and ANPCyT, Argentina. 
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